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Abstract 



We have modified the colour screening theory of Chu and Matsui by properly incorpo- 
rating bag model equation of state for quark gluon plasma (QGP). We have also chosen the 
pressure parametrization rather than parametrizing energy density in the transverse plane. 
We assume that the QGP dense medium is expanding in the longitudinal direction obeying 
Bjorken boost invariant scaling law. Sequential melting of ^ and J/V' is also considered 
in this scenario. We have applied above formulation to the recent PHENIX experimental 
data of J/tl' suppression in Au+Au collisions at RHIC. We find that the model gives a 
good description of data at mid-rapidity in terms of survival probability versus number 
of participants without any necessity of implementing (3+l)-dimensional expansion of the 
deconfined medium. 
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1 Introduction 



One of the main important aims of the Heavy-ion cohision experiments is to detect and find 
the properties of the state of matter in which quarks and gluons are deconfined as predicted 
by Lattice Quantum chromodynamics (QCD) above a critical temperature of the order of 
Tc ~ 0.17 GeV for a baryonic chemical potential hb = [Ij. The suppression of heavy 
quarkonia {J/ip, Xc ip , TT) due to colour screening analogous to Debye charge screening in 
QED plasma, has long been proposed as a probe of deconfinement in the dense partonic 
medium. Heavy quarkonia are thus considered as one of the most promising candidates to 
study the formation and the properties of QGP. In the deconfined state, the interaction 
between heavy quarks and antiquarks gets reduced due to colour Debye screening effects 
leading to a suppression in J/ip yields. Matsui and Satz [2j predict that the binding energy 
of the cc pair into J/V' mesons will be screened in the presence of a QGP medium, leading 
to the so called J/ip suppression. In the relativistic heavy ion collisions the J/'0 suppression 
has been recognized as an important tool to identify possible phase transition from Hadronic 
matter to quark-gluon plasma (QGP). The J/V' resonance states are produced at the initial 
(prethermal) stage of heavy ion collisions because of their large masses. Their small widths 
also make them insensitive to final state interactions. Therefore, their evolution probes the 
deconfined state of matter in the early stage of collisions [3]. 

Recently, high statistics data of Au+Au collisions at center-of-mass energy ^/sJlN = 200 
GeV at RHIC in Brookhaven National Laboratory have become available [4j. It is observed 
that J/V' yield in central Au+Au collisions is suppressed by a factor of nearly 4 at mid- 
rapidity and 5 at forward rapidity relative to that in p+p collisions scaled by the average 
number of binary collisions. Cold nuclear matter (CNM) effects such as nuclear absorption, 
shadowing and anti-shadowing are also expected to modify the J/ip yield. CNM effects due 
to the gluon shadowing and nuclear absorption of J/tp at the RHIC energy were evaluated 
from the J/ip measurement in d+Au collisions at RHIC [5j. PHENIX d+Au data show 
that CNM effects are smaller at RHIC than at SPS and can be reproduced by considering 
nuclear absorption cross-section of J/V' with nucleons and Au nucleus is of the order of 
3 mb and also incorporating a nuclear-shadowing effect which considers the depletion of low 
momentum partons in a nucleon embedded in a nucleus compared to their population in 
a free nucleon. Several groups have measured the J/^jJ yield in heavy- ion collisions (for a 
review of the data and interpretations see the refs. [HIT]). A J/'ip suppression obtained by 
the NA50 Collaboration at SPS [HUH [10] could be reproduced by various theoretical models. 
A larger suppression is expected to occur at RHIC compared to SPS due to a larger energy 
density present in the medium [111 I12| . However, the level of suppression is not much 
different from that observed by NA50 experiments. Model calculations assuming colour 
screening of J/V' state in a medium depict a much larger suppression at RHIC energies due 
to the presence of a large parton density, and a higher temperature as well as a large lifetime 
of the system. On the other hand, several models predict that J/ip yield will result from a 
balance between dissociation [131 [14] due to thermal gluons along with colour screening [15] 
and enhancement due to coalescence of uncorrelated cc pair |11[ [TBI [T7] which are produced 
abundantly in the initial stage of collisions at the RHIC energy |18[ [19] . 

We show in the present work that centrality (i.e., impact parameter or number of par- 
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ticipant nucleons Npart) dependence of the J/ip suppression in Au+Au collisions data at 
mid-rapidity recently observed by PHENIX experiment at RHIC can be explained invok- 
ing a QGP scenario only based on the bag model equation of state (EOS). We demonstrate 
that the survival probability pattern is well reproduced within present J/-^-!- hydro model by 
including the sequential melting of Xo "0 aiid J/ip in the longitudinally expanding plasma. 



2 Formulation 



Following closely the basic theme of refs. [15^ [20] but also highlighting the important dif- 
ferences at appropriate places, our formulation proceeds through the following five stages: 



2.1 Description of the medium 

In a plasma composed of u, d quarks and gluons, the temperature T(x), energy density 
e{x) and pressure P{x) depend on the time-space point x = {t, x) in the fireball frame. 
Assuming massless partons and local thermodynamic equilibrium, the bag model equation 
of state [21] reads 

e = aT^/ci + B ; P = aT^ - B ^ a{T^ - T^) (1) 

c„ = -:— ; a = ; B = ^ a 1^ 

de ' 90 ' 45 

where ~ 0.17 GeV is the critical temperature for hadron-QGP phase transition, ~ 1/3 
is the square of the velocity of sound in the medium, the coefficient a contains information 
about the degrees of freedom, and B is the bag constant. It should be emphasized that 
above choice of B gives a value B = 0.405 GeV/fm^. At T = Tc pressure of the QGP 
medium is zero and this means that QGP behaves as an ideal gas above this temperature. 
It would also be worthwhile to point out that Chu and Matsui ^IS] used bag model equation 
of state in their colour screening scenario to estimate the proper screening energy density 
yet they droped B while defining as the ratio P/e. For consistency we shall retain B 
throughout our analysis in this paper. 

The recent lattice QCD simulations do not support well the first order phase transition 
between the QGP and the hadron gas as is the case we get in the bag model EOS. Still, 
we want to emphasize that the bag model properly parametrizes many features of the EOS 
with a rapid change in the entropy density as a function of temperature and hence often 
continues to be an useful tool employed in the hydro dynamical calculations. If a QGP was 
produced in central collisions of identical nuclei and is expanding longitudinally, then local 
thermodynamic observables become function of the lateral coordinate r together with the 
proper time variable r defined by r = (t^ — z^)^/"^ . For algebraic convenience we define the 
dimensionless ratios 

T . T{T,r) e{T,r)-B ^_ P{T,r) + B 

Ti ' T{Ti,r) ' e{Ti,r)-B ' P{Ti,r) + B ^' 
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where Tj ~ 1 fm/c is the proper time for initial thermahzation of the fireball. Then Bjorken 
differential equation de/dr = — (e + P)/t leads to the scaling solutions or cooling laws 

l = p = f^ = f-i . q=l + cl (3) 

The cooling laws written by [151 120] do not exhibit the effect of B at all. 



2.2 Pressure profile 



Examination of (1) reveals that the pressure almost vanishes at the transition point Tc, 
i.e., becomes very small in the hadronic sector. This is also supported by the phase diagram 
given by Blaizot [22] . Hence on the transverse plane z = of the fireball and with ti = Ti 
we choose the pressure profile 

/ 2\l^ 

P{ti,r) = P{ti,0)h{r) ■ h{r)^\l-^\ e{RT-r) (4) 

where the coefficient P(ti, 0) is yet to be specified, Rt denotes the radius of the cylinder and 
it is related with the transverse overlap area Aj- of the colliding nuclei by Rt = \J AtIt^- 
The power /3 ~ 1 depends on the energy-deposition mechanism, and Q is the unit step 
function. Clearly, our pressure is maximum at the center of the plasma but vanishes at the 
edge Rt where hadronization occurs. In contrast, [151 120] used a similar parametrization 
for the energy density, ignoring the fact that e should suffer a jump by qB /c^^ ~ 45 across 
the phase transition point. The factor P(ti,0) is related to the mean pressure < P >i over 
the cross-section and to the corresponding average energy density < e >j via 

P{U,0) = {l+f3)<P >i= (1 + f3){cl <e>i-q B} (5) 



Both experimentally and theoretically the determination of J/V' survival probability 
S is of paramount importance. This quantity S, in principle, is a function of transverse 
momentum pT, rapidity y, and the impact parameter b (or number of participants Npart-) 
Assuming y ^ and 6 0, the original Chu and Matsui model [15] was proposed to explain 
the PT dependence of the J ftp suppression in terms of the time dilation of the formation 
time. However, in the present work our aim is to adopt the dilated formation time concept 
for describing centrality dependence of the J/^p suppression data at mid-rapidity which 
recently became available from PHENIX experiment at RHIC. In order to illustrate this, 
we take the initial average energy density < e >i in terms of the number of participating 
nucleons Npart [23j (which depends on the impact parameter 6), given by the modified 
Bjorken formula: 

<6>,= -^^ ; AT = 7rR'^ (6) 
Atti dy 

where At is the transverse overlap area of the colliding nuclei and cIEt / dy is the transverse 
energy deposited per unit rapidity of output hadrons. Both depend on the number of 
participants Npart [24J and thus provide centrality dependent initial average energy density 
< e >j in the transverse plane. Next, ^ is a phenomenological scaling factor discussed latter 
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in Sec. 3 in conjunction with the self-screened parton cascade model. Of course, [151 [20] do 
not write any factor in (6) so that their initial average energy density remains substantially 
underestimated in numerical applications so that they are also compelled to employ values 
determined by the self-screened parton cascade model. 



2.3 Constant pressure contour 



It is well known that a cc bound state kept in a thermal medium feels a colour screened 
Yukawa potential and it melts at the Debye temperature Tjj which corresponds to the 
energy density and pressure Pg given by 

Td>T, ; e, = ar^/c2 + S ; P, = aT^ - B (7) 

For any chosen instant t and on the z = plane the contour of constant pressure Pg is 
obtained by combining the cooling laws (3) with the profile shape (4) to yield 

P = Ps + B _ ~g 
~ Piti,0)h{r)+B 

Setting r = the maximum allowed tilde time tso (during which pressure drops to Pg at 
the center) can be identified as 



fP{ti,0) + B 



1/9 



with P{ti,0) read-off from (5). Thereby the said locus takes the more convenient form 

.2 \/3 



.X fi -B/iPs + B) , , 

1 t\ = Hs(t) = ^ ,\ \ (10) 

R't t%-B/{Ps + B) ^ > 



Our above result generalizes a similar expression derived by |15t |20] for the B = case. 



2.4 J/ip kinemitics and screening radius 



Consider an interacting cc created at i = at the location (r^, (p^) on the z = plane 

having mass m^, momentum p^, energy = ^ + , velocity = p^/p^, and dilation 

factor 7^ = p^/m^. In the fireball frame the pair will convert itself into the physical J/ip 
resonance after the lapse of time tp = j^Tp (with tf being the intrinsic formation time) 
provided the temperature T < Td- Prom the locus (10) we deduce the so called screening 
radius 

r, = Rt{1- H'/^{tF)}'/^9{l - Hsitp)} (11) 

which marks the boundary of the circular region where the quarkonium formation is pro- 
hibited. Hence, the pair will escape the deadly region and form quarkonium if 

\r^ + v^tF\>rs (12) 
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Analysis of the kinematic condition is greatly simplified if the J /ip is moving with pure 
transverse momentum pj- in the direction (i.e., mid-rapidity domain). Then, for escape the 
trigonometric condition to be obeyed becomes 

cos{^)>Y ■ [{rl-r^)m-rlpl/m] 
- 2rTFPT 

where all J/tp suffixes have been omitted without any loss of generality. 



2.5 Survival probability 



Suppose the radial probability distribution for the production of cc pair in hard collisions 
at r is ^ 

f{r) « - ^(^T -r) ; a ~ 1 (14) 

Then, in the colour screening scenario, the net survival probability for the quarkonium 
becomes 

where the maximum positive angle (pmax allowed by (13) is read-off from 

vr ify<-l 
. 7r-cos-i|y| ifo>y>-i 

' COS \y \ iio<y<i 
if y > 1 



Although the formulae (12-16) exist in the literature, yet we have found and corrected 
a serious misprint occurring in I20j concerning (pmax in the range — 1 < 1" < 0. In actual 
practice, it has been found that only about 60% of the observed J/ip originate directly in 
hard collisions while 30% of them come from the decay of Xc and 10% from the ip . Hence, 
the weighted survival probability of J /ip becomes 

5(pt) = 0.6 5^ + 0.3 5;,, + 0.15^' (16) 



3 Numerical Work 



Table 1 collects the values of various parameters used in our theory and the following expla- 
nations are relevant. The value Tc = 0.17 GeV is in accord with lattice QCD results. The 
choice = 1/3 is most common for free massless partons although for partons which carry 
thermal mass or interact among themselves may be different like 1/5 |20j . The selection 
P = 1 indicates that the energy deposited in the collision is proportional to the number of 
nucleon-nucleon encounters, i.e., to the nuclear thickness. Also, relevant properties of the 
various quarkonium species in a thermal medium are displayed in Table 2. It is clear that 
the dissociation temperature To gradually decreases in going from J/ip to Xc ip . 
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Table 1: Various parameters used in the theory. 

T, (GeV) cj B (GeV/fm^) (GeV/fm^) /? ^7" 
0.17 1/3 0.405 1 0.5 



Our numerical procedure proceeds through the fohowing steps: 

(i) Before finding the centrality (or impact parameter) dependence of J/ij: suppression it 
is necessary to know the initial average energy density < e >j in terms of the number of 
participants Npart- For this purpose, we extract the transverse overlap area At and the 
pseudo-rapidity distribution dET/dr] reported in ref. |24j at various values of number of 
participants Npart- These dET/dr] numbers are then multiplied by a constant Jacobian 1.25 
to yield the rapidity distribution dET/dy occurring in (6). 

(ii) The original Bjorken formula although provides an estimate of the initial energy density 
qualitatively but, unfortunately, it under-estimates the initial energy density which can 
cause the suppression of only Xc and ifj but not of J /il^- Hence, a scaling-up factor ^ = 5 
has been introduced in (6) in order to obtain the desired < e >i= 45 GeV/fm^ |25j for 
most central collision. The appropriate characterization of kinematic quantities in Au+Au 
collisions is presented in Table 3. The relatively large values of our < e >j have the following 
justification: 

These are consistent with the predictions of the self-screened parton cascade model [26], 
these agree with the requirements of hydrodynamic simulation [25] which fit the pseudo- 
rapidity distribution of charged particle multiplicity dNch / drj for various centralities already 
observed at RHIC, and these can cause melting of all the quarkonium species listed in Table 
2. 

(iii) Next, we calculate the time tso for the pressure to drop to P, at the origin and thereby 
deduce the screening radius with the help of (9, 10). 

(iv) Next, the quantity Y is computed from (13) which sets the condition for the quarkonium 
to escape from the screening region, and the limiting values of the (pmaxii^) are constructed 
using equation written just below (15). 

(v) Finally, the survival probability S{pt), at specified pT but varying Npart is evaluated 
by Simpson quadrature from (15,16). 

(vi) In order to compare the above analysis with the actual experiments it is necessary to 
convert the J/V' suppression data available in terms of nuclear modification factor Raa 
into the accepted def. of survival probability S{pt) [271 ESI 129] namely 

where Raa is the standard nuclear modification factor and Raa^ i^ ^ contribution to Raa 
originating from CNM effects constrained by the data of d -|- Au collisions. 

Now we turn to physical interpretations of our results. 
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S(p^) of J/psi including sequential decay at RHIC 



0.8- 



^ Present Model 

• PHENIX Au+Au (|y|<0.35) 

^ Chu and Matsui Model 
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Figure 1: Graph depicting the survival probabiHty versus number of participants at fixed 
PT = 2 GeV/c ahowed by invariant pr spectrum of J/i/j measured by PHENIX H EH [23 
when sequential melting is considered. 
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S(p.^) of directly produced J/psi at RHIC 
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Figure 2: Same as Fig. 1 except that sequential melting is not considered. 



8 



Table 2: Masses, formation times and dissociation temperatures of J/i)-, Xc V'' [30j . 





i, 


Xc 


^' 


m (GeV) 


3.1 


3.5 


3.7 


Tp (fm) 


0.89 


2.0 


1.5 




2.1 


1.16 


1.12 



Table 3: Kinematic characterization of Au+Au collisions at RHIC [3] 



< e >i (GeV/fm^) i?^ (fm) 



Nuclei A/sjVjv (GeV) g 



part 



Au+Au 200 



5.0 



22.0 


5.86 


3.45 


30.2 


7.92 


3.61 


40.2 


10.14 


3.79 


52.5 


12.76 


3.96 


66.7 


15.69 


4.16 


83.3 


18.58 


4.37 


103.0 


21.36 


4.61 


125.0 


24.38 


4.85 


151.0 


27.37 


5.12 


181.0 


30.52 


5.38 


215.0 


34.17 


5.64 


254.0 


37.39 


5.97 


300.0 


41.08 


6.31 


353.0 


45.09 


6.68 
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4 Results and Discussions 



We present our numerical results under the following two headings: 

Sequential decay included. Figure 1 shows the variation of Survival probability S{pt) with 
respect to number of participant nucleons Npart at transverse momentum px = 2 GeV/c 
allowed by the invariant pt spectrum of J /ip measured by PHENIX [3] at RHIC energy by 
including sequential decay of Xc and ip into Jip calculated from (15, 16). The three curves 
correspond to experimental data [3; solid circles], our formulation using the bag model EOS 
(solid triangles) and old Chu-Matsui theory without bag model EOS [18; solid rhombus]. 
It is obvious from Figure 1 that S{pt) decreases with increase in number of participant 
nucleons in all the three curves because of the growth of the initial energy density. However, 
the agreement between experiment and our formulation is very good (characterized by the 
gradual fall of S{pt) with Npart) which justifies our use of the pressure parametrization 
containing bag constant B. On the contrary, the agreement between experiment and Chu- 
Matsui theory is poor (as characterized by sharp fall of S{pt) with Npart) which negates 
their use of the energy parametrization ignoring B. 

An important remark must be added at this juncture. Although the quality of fit between 
our model and experimental data is excellent for central collisions (i.e., relatively large values 
of the Npart) it worsens somewhat for non-central collisions (i.e., smaller values of Npart)- 
This may be due to the fact that our model is based on the assumption of Bjorken boost 
invariant longitudinal expansion, which requires cylindrical symmetry about the collision 
axis. This cylindrical symmetry may be violated up to some extent for most non-central 
collisions. 

Sequential decay excluded. For the sake of comparison Figure 2 shows similar plots when 
the theory considers only the suppression of directly produced J /ip^s. The old Chu-Matsui 
model again shows much larger J /tp suppression (lower values of S{pt)) for central as well 
as non-central collisions as compared to experimental results. In contrast, the agreement 
between the trends (i.e., shapes) of our theoretical curve and experimental data is much 
better although their absolute magnitudes differ noticeably for central as well as non-central 
collisions. Thus, inclusion of the sequential melting of the charmonia ^23] seems to be 
essential for fitting the data. 

5 Summary and Conclusions 

We have analyzed centrality dependence of the J/^l^ suppression data in Au+Au collisions 
available in terms of survival probability versus number of participants at mid-rapidity 
from PHENIX experiment at RHIC [31 EH [28l EHj by using the modified Chu and Mat- 
sui [15] model with J/^/^-|- hydro framework. We have simultaneously incorporated several 
features namely, bag model equation of state, longitudinal expansion of the QGP, pressure 
parametrization (rather than energy parametrization) in the transverse plane, scaling factor 
in Bjorken formula to generate higher energy densities (compatible with hydro dynamical 
simulations) , dilated J/ ip formation time and sequential melting of charmonia in the colour 
screening scenario. 
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In conclusion, our present JZ-i/^+hydro formulation presents a very good agreement 
with the experimental mid-rapidity data of J /tp suppression versus number of participants 
observed by PHENIX experiment at RHIC. Of course, a more rigorous treatment of the 
research problem would require incorporation of additional complications such as (3+1)- 
dimensional expansion of the QGP [20^ I27j. EOS built from non-equilibrium fugacities of 
the partons [31], gluonic dissociation [13^ I14j. etc. Work on these topics will be taken up in 
a future communication. 
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